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Abstract 

The symmetrized density matrix renormalization group approach is applied 
within the extended Hubbard- Peierls model (with parameters U/t, V/t, and bond 
alternation 5) to study the ordering of the lowest one-photon (l 1 !?") and two-photon 
(2 1 A+) states in one- dimensional conjugated systems with chain lengths, N, up to 
N = 80 sites. Three different types of crossovers are studied, as a function of U/t, 
8, and N. The [/-crossover emphasizes the larger ionic character of the 2A g state 
compared to the lowest triplet excitation. The 5 crossover shows strong dependence 
on both N and U/t. The ^-crossover illustrates the more localized nature of the 
2A g excitation relative to the 1B U excitation at intermediate correlation strengths. 

PACS Numbers: 71.10.-w, 71.24.+q, 78.66. Qn 

Recently, much attention has focused on the luminescence properties of conjugated 
organic materials because of their potential for application in display devices [1]. These 
studies have underscored the importance of the structure of low-lying electronic excited 
states. Specifically, a major parameter is the relative ordering of the lowest dipole allowed 
singlet (l 1 !?") state and the lowest dipole forbidden singlet (2 1 A+) state, in the light of 
Kasha's rule which relates molecular fluorescence to the lowest excited singlet state. 

It is well established that correlated electron systems behave differently from indepen- 
dent electron systems, especially in the case of excitations. Earlier work has shown that 
the lowest optically forbidden excited state 2A g lies below the optically allowed excited 
state 1B U in polyene molecules [2] ( thus preventing any significant luminescence in such 
compounds ), while an independent electron model gives the opposite picture; similar 
results have been found by Periasamy et al. in the case of polycrystalline sexithienyl [3] 
or Lawrence et al. in single crystal polydiacetylene [4]. These examples serve as evident 
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manifestation of electron correlation in conjugated molecules. The influence of electron 
correlation has also been considered as the main origin of lattice dimerization leading to 
the view that conjugated polymers are rather Mott insulators than Peierls insulators [5]. 
In the context of third-order nonlinear optical response and photo-induced absorption, 
the role of higher-lying excited states derived from correlated electron models has also 
been emphasized [6]. 

It is also important to stress that the electronic and optical properties of conjugated 
oligomers and polymers differ depending on whether the compounds are in the gas phase, 
in solution, or in the solid state [7]. The chemical environment affects the geometric 
structure as well as the electron correlation strength, the latter via dielectric screening [8] . 
Furthermore, the characteristics of the conjugation defects present in oligomers depend 
on chain length, which emphasizes the influence of quantum size effects [9]. 

In view of these features, we believe that these three factors: (i) geometric structure, 
(ii) strength of electron correlation, and (iii) quantum confinement, are most relevant for 
the study of the photo- and electro-luminescence response in organic conjugated chains. 

Previous studies of the 1B/2A crossover behavior have been carried out for short chain 
systems, so that in the independent electron limit, the 2A energy is significantly higher 
than that of IB due to the discreteness of the molecular orbital energy spectrum. In this 
zero U limit, the 2A state corresponds to single HOMO to LUMO — 1 (or HOMO + 1 to 
LUMO) excitation while IB is a HOMO to LUMO excitation. According to previous 
results [10], as electron correlation U is turned on, the gap between the ground state and 
the 2A state narrows while the gap to the IB state increases; the states thus cross at a 
given Hubbard correlation strength U c . This, we refer to as the [/-crossover. However, 
for an infinite chain, the 2A and IB states both occur at the same energy in the H'uckel 
limit (U = 0). If the 2 A and IB states were evolving in a manner identical to that in the 
short chains, these states would never cross with increasing U. Thus, for a given U, their 
must occur a crossover from the short chain behavior to the long chain behavior; this we 
refer to as the iV-crossover. 

It was noted by Soos, Ramasesha and Galv'ao [11] from exact diagonalization studies 
of short chains that a similar crossover occurred with variation of the bond-alternation 
parameter 5, which we refer to as the 5- crossover. The 5-crossover was studied by moni- 
toring the optical gap and the lowest singlet-triplet ( spin ) gap; the critical 5 C for a given 



2 



correlation strength was determined by the value of 5 at which the optical gap equals 
twice the spin gap. These authors further described the system as behaving band-like for 
5 values above 5 C and correlated-like for 5 values below 5 C . However, as was pointed out 
in Ref. [12], increasing bond alternation does not lead to the band picture, because the 
binding energy of the IB exciton increases with increasing 5, an obvious indication that 
electron correlation increases at the same time. 

In this Letter, we present a thorough study that encompasses the three kinds of 
crossovers, namely the U, N, and 5 crossovers in conjugated chains, by employing the 
symmetrized density matrix renormalization group ( SDMRG ) theory. The SDMRG ap- 
proach [13] is currently the most reliable many-body method for calculating the low-lying 
excited states with high accuracy for relatively large systems and for a wide range of 
model parameters. The model Hamiltonian in this study is the extended Hubbard- Peierls 
Hamiltonian which reads: 

<ij>,a i <ij> 

where rji — [1 — S(— 1)*], 5 is the dimensionless dimerization parameter, U is the on- 
site Hubbard repulsion ( in units of t, the nearest-neighbor hopping integral ) and V is 
the nearest-neighbor charge density-charge density interaction. The 5 term serves as a 
structural parameter in the simplest way, if we assume linear electron-lattice coupling in 
the static limit; as has been pointed out before, the V-term is crucial to the understanding 
of the optical excitation spectrum, namely the excitonic effect [14]. The present model 
can be regarded as the minimal correlated model for conjugated systems. Note that, the 
meaningful phase corresponds to the BOW ( bond-order wave ) regime, namely, V < U/2 
[15]. By comparing experimental data for a series of polyene molecules, we find that the 
parameter set (t,U,V,5) = ( 2.4eV, 7.2eV= 3* , 0.4 U, 0.07) gives the best fit for the IB, 
2 A, and even higher energy A (mA) states [16]. We thus set V/U — 0.4 without losing 
generality. 

The density matrix renormalization group method is the most accurate numerical 
method for determining the ground and low-lying excited states of quasi-one-dimensional 
correlated electron systems with short range interactions [17]. In the usual DMRG proce- 
dure, it is difficult to target the 1B U state as there are many states that appear between 
it and the ground state, with the number of these states increasing with U and chain 
length N. However, in a symmetrized DMRG technique that exploits spin parity, C 2 
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symmetry, and electron- hole symmetry, the 1B U state is the lowest state in the subspace 
e B~ . Incorporating these three symmetries thus allows us to determine the 1B~ and the 
2A/ j state energies with unprecedented accuracy for chains of up to 80 sites. We choose 
to truncate the space of density matrix eigenstates to 100 ( m = 100 ) in most cases. 
For smaller U and 5, however, we choose a larger value of m (= 150) in order to achieve 
consistent accuracy. 

We contrast the U — crossover for short ( N — 8 ) and long ( N = 80 ) chains for fixed 
alternation 5 = 0.07 in Fig. I. It is well known that in the strong correlation limit, the 
2A state becomes a spin excitation which is gapless in the limit S = and this state can 
be described as composed of two triplets. Thus, increase in correlation strength should 
show a decrease in the 2A energy [18]. However, we note that in the N = 8 chain, the 
two-photon state energy remains nearly constant before decreasing for values of U / 1 larger 
than 2.0. In the longer chain, the 2A g energy increases even more rapidly with increasing 
correlation strength, than the 1B U energy. This implies a substantial ionic contribution to 
the 2A g state in long chains besides the covalent triplet-triplet contribution. This result 
constitutes the first clear illustration of the importance of quantum size effects. We find, 
however, that the critical correlation strength, U c , at which the crossover occurs is nearly 
independent of the chain length N; in both N = 8 and iV = 80 cases, U c is around 2.5t. 

For fixed correlation strength ( U/t — 3 and 4 ), we present the (^-crossover results 
for N = 8 and 80 in Fig. 2. We find that the critical 5 value, 5 C , strongly depends on 
chain length. For U/t = 3, the S c values are found to be 0.15 and 0.09 for N = 8 and 
80, respectively; for U/t = 4, they are 0.32 and 0.22. Thus, S c has both strong iV and 
U dependence. We also show in Fig. 2 the crossover behavior between the 1B U energy 
and twice the lowest triplet energy, E T . This crossover occurs at systematically smaller 
5 values, again emphasizing the larger ionic character present in the 2A g state compared 
to the lowest triplet state. 

Most interestingly, we find one more crossover behavior, which is the iV-crossover, in 
the case of intermediate U/t and medium to large 5 values. We observe that the 1B U and 
2A g states cross over for fixed U/t and S as a function of N, the chain length. The critical 
lengths are actually fairly insensitive to U and 8. In Figs. 3a ( U/t — 3, 5 — 0.12 ) and 3b 
{U/t = 4.0, 5 = 0.27 ), we find this crossover for iV = 14 and iV = 12, respectively. This 
is a direct theoretical observation of quantum confinement induced crossover. It is related 
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to the fact that the 2A g excitation is more local in character with a shorter characteristic 
length than the 1B U state. Thus, the \B U excitation is stabilized over longer length 
scales than the 2A g excitation. This is seen as a more rapid saturation in the 2A g energy- 
compared to the 1B U energy, as a function of chain length. We note that this crossover 
can also be seen from Fig. 2 where the S c values show a decrease in going from N = 8 
to N = 80. This behavior can only exist for intermediate correlation strength: for weak 
correlation, there does not exist any crossover and 2A g lies above the 1B U state for all 
chain lengths as seen from Fig.l; at large values of U/t, we are in the atomic limit, a 
crossover is not expected and the quantum size effects are largely suppressed. It has been 
widely accepted that the conjugated molecules fall in the intermediate correlation regime; 
thus, the confinement-induced crossover is realistic. 

It is most relevant in this context to stress that, for smaller 5 values, even though 
no crossover occurs, the 1B/2A gap decreases in long chains. When considering the 
parameters adapted to polyenes ( U/t = 3, V = OAU, 5 = 0.07 ), both 1B U and 2A g state 
energies decrease with N, but 2A g state saturates in a way much faster than IB U , because 
the former ( mostly covalent ) is more localized than the latter. Consequently, it would 
be inappropriate to extrapolate the 1B/2A gap in polyacetylene from data obtained on 
short chains [19]; our results imply that the actual gap is significantly smaller than the 
result of such an extrapolation. 

To conclude, we have employed the accurate numerical density matrix renormaliza- 
tion group technique with symmetry adaptation to study the ordering of the lowest one- 
photon and two-photon states in conjugated oligomers and polymers within an extended 
Hubbard- Peierls model. Three kinds of crossover, namely a [/-crossover, a (^-crossover, 
and a iV-crossover, have been demonstrated. The iV-crossover is related to quantum fi- 
nite size effects and crucially depends on the characteristic length of the excitations; this 
characteristic length is of the same order as the chain length at intermediate correlation 
strengths and degrees of dimerization. 
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Figure Captions 

Figure 1. Crossover on U for 5 = 0.07. 

Figure 2. Crossover on 5 for (a) U/t = 3; (b) U/t = 4. 

Figure 3. Crossover on N for (a) U/t = 3 and 5 = 0.12; (b) U/t = 4 and 5 = 0.27. 
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